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ABSTOACr 


The  models  for  determining  convergence  zone  Gain  (G)  were  developed  using 
a  linearized  Sound  Speed  Profile  (SSP)  and  applying  ray  tracing  theory.  The  SSP 
was  divided  into  three  cases;  bilinear,  bilinear  with  isospeed  layer,  and  bilinear  with 
mixed  layer.  Two  analytical  solutions  were  developed  using  Taylor  series  and 
binomial  series  expansions  to  determine  G,  one  for  the  bilinear  and  bilinear  with 
isospeed  layer,  and  the  other  for  the  bilinear  with  a  mixed  layer.  The  solutions  for 
G  are  exclusively  a  function  of  the  SSP  gradients.  Each  solution  was  compared  to 
the  solutions  from  ray  tracing  and  the  solutions  from  the  Integrated  Carrier  ASW 
Prediction  System  (ICAPS)  (which  runs  on  mainframe  computers  and  requires  more 
data  in  addition  to  the  SSP).  When  the  SSFs  were  not  too  unusual,  the  solutions  for 
G  were  fairly  close  when  compared  to  ray  tracing  and  ICAPS. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Convergence  .zone  (CZ)  parameters  such  as  range  to  CZ,  CZ  width,  and  CZ 
gain  are  typically  obtained  using  large  mainframe  computers  implementing  complex 
software  programs.  The  U.S.  Navy  utilizes  the  r’leet  Numerical  Oceanographic 
Center  (FNOC),  Monterey,  California,  to  predict  these  parameters.  The  frequencies 
of  interest,  source  and  receiver  depths,  Sound  Speed  Profile  (SSP),  wave  height,  time 
of  day,  date,  and  geographic  location  are  used  in  obtaining  the  CZ  parameters. 

In  addition  to  FNOC,  the  fleet  may  use  the  Integrated  Carrier  ASW  Prediction 
System  (ICAPS).  ICAPS  was  designed  for  use  onboard  ships.  It  consists  of  several 
databases  of  environmental  data  and  can  provide  the  user  with  a  cornucopia  of 
acoustic  parameters. 

Both  of  these  methods  require  the  use  of  mainframe  or  "mainframe-like" 
computers.  Several  programs  for  hand-held  and  desktop  calculators  have  been 
written  in  attempts  to  obtain  quick  estimates  of  CZ  parameters  based  on  current 
bathythermograph  (BT)  data.^  However,  none  of  the  calculator  programs  is  usable 
on  IBM-based  disk  operating  system  (MS-DOS)  running  on  personal  computers 
(PCs).  With  the  advent  of  notebook  PC’s,  and  with  the  palmtop  PC  on  the  horizon, 
a  definite  need  exists  for  models  which  will  run  on  an  MS-DOS  based  computer. 

^  See  the  Bibliography  for  further  information. 
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B.  DESIRABLE  QUALITIES  FOR  COMPUTER  MODELS 

The  following  is  a  list  of  general  guidelines  for  computer  models 


[Ref.  I:pp.l3,14]. 

•  A  minimum  amount  of  data  should  be  required  to  run  the  program. 
Additionally,  any  data  used  should  not  be  entered  more  than  once. 

•  Output  should  give  range  to  the  inner  edge  of  the  CZ,  width  of  the  CZ,  and 
CZ  gain.  The  output  should  also  give  a  repeat  of  the  inputs. 

•  The  program  should  run  quickly. 

•  The  program  should  be  based  on  acoustic  theory  rather  than  general  thumb 
rules  and  agree  with  mainframe  prediction  models. 

C.  PURPOSE 

The  purpose  of  this  study  is  to  provide  a  MS-DOS  based  program  following  the 
general  guideline  for  desirable  qualities  for  PC  models  listed  above.  The  model  will 
be  based  on  ray  tracing  vice  ray  mode  theory. 
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n.  CZ  RAY  TRACING  THEORY 


A.  WHY  RAY  TRACING? 

To  simplify  the  mathematics  and  programming  involved,  ray  tracing  theory  was 
selected  over  normal  mode  theory.  Ray  tracing  provides  an  easy  visualization  of  how 
sound  energy  propagates  through  the  ocean,  even  with  complicated  SSP  and  bottom 
structures.  A  sound  ray  can  be  traced  from  the  source  to  the  desired  range  without 
any  need  for  complex  calculations.  As  diffraction  of  the  sound  energy  becomes  more 
noticeable  at  lower  frequencies,  and  as  ray  tracing  theory  does  not  consider 
diffraction,  ray  tracing  becomes  more  valid  as  the  frequency  increases. 

Normal  mode  theory  gives  a  complete  description  of  the  sound  field.  However, 
it  cannot  easily  handle  real  boundary  conditions.  All  frequencies  can  be  used  but 
complexity  increases  with  the  frequency.  [Ref.  2:p.l22] 

B.  CX)NVERGENCE  ZONES 

Imagine  the  sound  rays  projected  from  a  source  at  the  surface.  Some  of  these 
rays  may  get  trapped  in  a  surface  duct  or  interact  with  the  bottom.  Other  rays  may 
be  able  to  penetrate  below  the  SOFAR  axis,  get  refracted  upward,  and  reacquire  the 
surface.  Starting  from  the  horizontal  and  increasing  the  angle  of  depression,  the  first 
ray  which  crosses  the  sound  axis  and  reacquires  the  surface  does  so  at  a  range  equal 
to  router .  As  the  depression  angle  increases  further,  the  range  at  which  the  ray 
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resurfaces  first  decreases  to  a  minimum,  then  increases,  resweeping  back  beyond 
router  For  omnidirectional  source  the  reswept  area  forms  an  annulus  about  the 
source  and  is  acoustically  brighter  than  the  surrounding  areas.  When  this  brightened 
area  is  compared  to  a  reference  area  at  1  meter  from  the  source,  the  CZ  gain  can 
be  calculated. 

Figure  1  illustrates  the  horizontal  distance  a  ray  travels  in  a  linear  sound  speed 
gradient  g. 


Figure  1  -  Horizontal  Distances 

For  the  case  shown,  the  sound  speed,  c,  is  expressed  as  a  function  of  depth  by 
c(z)  =  Co(l+gz).  At  depth  Zj,  the  sound  speed  is  Cj,  and  the  downward  angle  of  the 
ray  is  Sj.  At  depth  z,,  the  sound  speed  is  c,,  and  the  downward  angle  of  the  ray  is 
02.  Snell’s  law  states  that  c/(g  cos0)  is  constant  along  the  ray  path. 

The  gradient  is  given  by 
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g= 


and  radius  of  curvature  R  of  the  ray  is 

R=_i_= — L_= — L_ 

gCOS0  gCOS0j  gCOS02 


The  horizontal  distance  Ar  the  ray  travels  and  the  change  in  depth  Az  are  given 
geometrically  by 


Ar=r,-rj  =R(sin02-sin0j) 


AZ=Z2-Zj=R(COS0j  -COS02) 

The  ordering  of  subscripts  in  the  above  equations  result  from  the  convention 
of  measuring  an  angle  of  elevation  positive  upwards  but  depth  is  defined  as  positive 
downward. 

1.  Range  to  Outer  Edge  of  the  Reswept  Zone 

The  ray  with  initial  depression  angle  of  0ou,  defines  the  range,  rou,„ ,  to  the 
outer  edge  of  the  reswept  zone.  At  this  range,  the  acoustic  energy  experiences  a 
sharp  drop  in  intensity.  For  the  bilinear  profile  illustrated  in  Fig.  2,  ©ou,  =  0°. 


5 


6 


Bottom 


2.  Range  to  Inner  Edge  of  the  Reswept  Zone 

The  ray  with  a  depression  angle  of  6in  defines  the  range,  rin^^  to  the  inner 
edge  of  the  reswept  zone.  Because  of  the  strong  convergence  of  neighboring  rays, 
with  the  formation  of  a  caustic  (where  sound  rays  cross  each  other)  at  the  inner  edge 
of  the  CZ,  there  is  a  strong  enhancement  of  the  intensity.  Additionally,  ray  theory 
does  not  provide  a  suitable  description  of  this  part  of  the  problem. 

3.  Range  to  the  CZ 

The  range  to  the  CZ  zone  will  be  defined  as  the  range  rou,„  [Ref.  1  :p.51]. 
For  the  remainder  of  the  discussion,  the  reswept  zone  will  be  considered  the  CZ,  and 
bottom  interactions  will  not  be  considered.  This  follows  a  consistent  course  as 
developed  in  earlier  models  which  have  also  ignored  bottom  interactions  [Ref.l;p.51]. 

4.  Width  of  the  Convergence  Zone 

The  CZ  width,  Ar^^  » is  the  distance  from  the  outer  edge  of  the  reswept 
zone  to  the  inner  edge  of  the  reswept  zone.  The  width  of  the  CZ  is  thus  taken  to 
be 

^''cZ"^ouicr“^inner 

5.  Convergence  Zone  Gain 

The  general  definition  for  Transmission  Loss  (TL)  is 


TL=101og 


fill)) 

U(r)J 
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where  1(1)  is  the  intensity  of  the  sound  beam  at  r  =  1  meter,  and  I(r)  is  the  intensity 
of  the  sound  beam  at  r  meters  from  the  source. 

By  conservation  of  energy,  the  following  relationships  for  the  intensities 
and  areas  through  which  the  rays  pass  normally  can  be  made. 

I(l)A(l)=I(r)A(r) 


where  A(l)  is  the  cross-sectional  surface  area  for  the  sound  beam  as  illustrated  in 
Fig.  3,  and  A(r)  is  the  cross-sectional  area  for  the  beam  at  r.  The  geometric  TL  is 
now 


TLg  =  101og 


(m] 


Figure  3  -  Reference  Area  Diagram  for 
CZ  Gain 
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The  area  A(l)  is 


A(l)=27r  A6cos0j„ 


Figure  4  pictures  a  cross  sectional  view  of  the  area  contained  by  the  CZ. 
The  intensity  will  be  assumed  to  be  uniform  over  this  area.  The  area  normal  to  the 
arriving  rays  of  the  CZ  annulus  [Ref.  l:p.62]  is 

^(^nncr)  “^^^inncr^ 


Figure  4  -  CZ  reference  area  (cross  section) 


9 


TLg  becomes, 


2ir  A0  cosej„ 


27rr:, 


inner 


h 


Ar 


cz 


=tane;„ 


Ar 


cz 


=sin8|„ 


Ar, 


h  =  _:^(l-cos0i„) 
sin0:„' 


When  the  small  angle  approximation  of  A 0  =  0i„  <<  1  is  substituted,  and 

A(l)=2:re„ 


TLg  is  now 
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TLg=201ogri„„„-Hl01og 

CZ  Gain  (G)  is  the  difference  between  the  TLg  from  spherical  spreading 
[201og(r)]  and  what  is  actually  observed  for  the  geometrical  spreading  loss.  (There 
are  two  definitions  for  G.  The  first  definition  includes  additional  losses  such  as 
absorption  in  G,  and  the  second  definition  considers  absorptive  losses  separately 
[Ref.  3:p,4C8]  [Ref.  4:p.7.1].  As  a  frequency  independent 
model  will  be  developed,  the  second  definition  for  G  will  be  used.) 

So  G  is  now 

G  =  101ogi!ill!li^  (2) 

In  real  ocean  conditions,  the  intensity  within  the  CZ  has  an  appearance 
similar  to  a  Rayleigh  distribution.  The  intensity  rises  sharply  from  zero  at  rj^ner  and 
gradually  declines  to  zero  at  rou,„  .  Because  the  exact  distribution  is  unknown  and 
a  single  number  is  desired  for  the  CZ  gain,  G  will  be  assumed  to  be  uniformly 
distributed  over  Ar^. ,  and  (2)  will  give  the  average  value  of  G.  Simplifying  G  to  an 
average  value  follows  earlier  work  done  in  similar  models  [Ref.l:p.40]. 
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m.  MODEL  DEVELOPMENT 


A.  GENERAL 

Several  approximations  are  useu  throughout  the  development  to  reduce 
equations  to  a  manageable  and  easily  understandable  form.  The  approximations 
used  are  contained  in  Appendix  A. 

Tlie  convention  used  when  labeling  angles,  sound  speeds,  and  radii  of  curvature 
is  the  subscript  is  always  tied  to  the  gradient.  As  illustrated  in  Fig.  5,  for  the  deep 
gradient  gj,  the  ray  entering  this  region  at  a  depth  Zj  has  a  radius  of  curvature  Rj  and 
an  initial  angle  entering  this  region  of  Bj. 

Sound 

Speed 


Figure  5  -  Reference  diagram  for  subscript  conventions. 
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The  three  SSP’s  which  will  be  examined  are  bilinear,  bilinear  with  isospeed 
layer,  and  bilinear  with  mixed  layer.  Figure  6  illustrates  the  three  profile  types. 


Figure  6  -  Bilinear,  Bilinear  with  Isospeed  Layer,  and  Bilinear  with  Mbced 
Layer  SSP’s 


B.  BILINEAR  SSP 

1.  Range  to  the  Outer  Edge  of  the  CZ 

The  ray  reacquiring  the  surface  at  rou,er  has  an  initial  depression  angle 

Qou,  =  0-. 

2.  Range  to  the  Inner  Edge  of  the  CZ 

Examination  of  Fig.  5  shows  that  the  range  from  the  source  (at  the 
surface)  to  where  a  ray  with  a  depression  angle  Bq  reacquires  the  surface  can  be 
expressed  as: 
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(3) 


r=2 


Ro(sin0j-sin0(j)  +  RjSin0j 


Using  Snell’s  Law: 


Co  Cj 


cos  00  COS0J 

The  radius  of  curvature  for  the  lower  layer  Rj  can  be  expressed  in  terms 
of  the  radius  of  curvature  for  the  upper  layer  Rq, 


and  when  Rj  is  eliminated  from  the  range  equation  (3), 


r=2R 


l+?2 


sin0j-sin0o 


Again  using  Snell’s  Law,  a  relationship  can  be  found  for  the  0o  and  0i , 


sin0j=^l-cos'0j 
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Manipulation  of  (the  above  equation)  to  obtain  rjnn„  is  cumbersome  but 
direct.  Differentiating  r  with  respect  to  6o  and  imposing 


dr 

dOo 


enables  0jn  to  be  solved  analytically  and  directly  input  into  (3)  to  solve  for  rjnnef . 

Without  further  simplification  the  solution  for  0o  yields  a  difficult 
transcendental  equation  for  0in .  However,  because  0o  is  generally  significantly  less 
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than  20*  (0.3491  radian)  the  trigonometric  functions  can  be  replaced  by  Taylor  series 
expansions  and  terms  of  third  and  higher  order  can  be  discarded.  By  replacing  the 
trigonometric  functions,  (4)  simplifies  to 


1 


l-ie; 


\2 


Further  manipulation  of  the  range  equation  using  the  binomial  expansions 
in  Appendbc  A  yields 


/'  'I 

2 

1  +  1° 

'>  Cl 

1  gj 

.2 

Co 

(5) 


Convenient  parameters  can  be  defined: 


(6) 


(7) 


and, 
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r=i+^ 

gi 


(8) 


The  final  form  of  the  range  equation 


r 


(9) 


Now,  differentiation  yields 


^100  go 


©or 


/I'Y'+eS 


Note  that  for  Bq  =  0,  dr/d0o  is  negative,  so  the  range  is  decreasing  with 
an  increasing  Bq  .  The  minimum  range,  ri„„„  is  found  by  determining  the  value  of 
00  =  Bin  for  which  dr/dBo  vanishes. 


dr 

=  2^® 

0:„r 

in 

dBo 

eo'0.n 

go 

=0 


Manipulating  this  equation  gives 
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and  finally, 


6fn(r'-l)-(Y'-l)'0 


1-Y- 

r’-i 


(10) 


Substituting  into  (9)  yields  an  finn^  of 


^inner 


go 


N 


a^  +  . 


r--i 


r-1 


Reducing  we  have 


^inner  ^ ' 


a" 


goN  r’--l 


!r=-l] 


and  finally, 
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(H) 

3.  CZ  Gain 

With  substitution  of  and  rou,„  into  (1),  the  CZ  width  can  be  expressed 


as 


Arcz=2J 

go 


r{^ -\/a-(r--l) 


Ar 


CZ 


=2 


go 


(12) 


Substitution  into  the  CZ  gain,  equation  (2),  gives 


G  =  101og 


Observe  that  G  is  only  a  function  of  the  upper  and  lower  gradients  and  the 
lack  of  dependance  upon  sound  speed,  depth  (although  these  are  Implicitly  contained 
in  the  gradients),  and  the  ray  departure  angles. 


4r 
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C.  BILINEAR  WITH  ISOSPEED  LAYER  SSP 

1.  Range  to  the  Outer  Edge  of  the  CZ 

As  with  the  previous  SSP,  r^ujcr  is  determined  by  using  an  initial  depression 
angle  of  6o  =  0°. 

2.  Range  to  the  Inner  Edge  of  the  CZ 

Adding  an  isospeed  layer  of  thickness  h  to  the  bilinear  SSP  gives  a  new 
range  equation 


r  =  2 


RQ(sin6j-sin0o)+hcot0j+RjSin0j 


(14) 


Manipulation  of  the  R’s  and  trigonometric  functions  yields 


r  Co 

r 

1  -y2  -tan0ft 

+  _ 

yh 

2  go 

n.  *  ^ 

COS-00 

1  2 

— 

N 

COS- 00 

Making  use  of  series  expansions  from  Appendbc  A, 


(15) 


An  additional  definition  which  will  prove  to  be  useful  is, 


(16) 


To  get 


dr 


eo-Oin 


=  0 


Which  gives  the  following  equation, 


0  = 


and  by  factoring  out  the  radicals, 


,  ,  .  ,  2erer„ 

a-*0r„=r=er„-_^" 


and  finally. 


o=e‘,[r=-i]+e;„Q=[2r^-3+2£r)*er,a- 


-  2er  e- 
r'-3+ — +— 

'>  4 

a 


(17) 
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Note  that  this  equation  is  cubic  in  Q-J  whereas  the  bilinear  case  the 
corresponding  equation  was  linear  in  6j„.  Additionally,  as  h  approaches  zero,  this 
reduces  to  equation  (10).  To  reduce  this  cubic  equation  to  a  more  manageable  form, 
the  relative  size  of  the  coefficients  for  will  be  considered.  Using  historical  data 
on  SSP’s  from  the  Atlantic,  Pacific,  and  Mediterranean  the  following  can  be  related 
[Ref.  l:pp.30,67]: 

Y<1 

a<l 

go<l 

0<h<250 
Co»1500 
Thus,  e<l 


Because  e  is  much  less  than  1,  the  e’  term  will  be  discarded.  By 
discarding  this  term,  an  equation  which  is  quadratic  with  respect  to  6,^  results  and 
can  be  expressed: 


o=e:„[r=-i]+era= 


r--2- 


2re 


-a’ 


(18) 


And  finally, 
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i2 


(19) 


1 
2 

— — —  -r--2-^+  r2-2-^  +4{r2-i) 

n|  2{r2-i)  [  i  0(2  J  a2  J 

3.  Width  of  CZ  and  CZ  Gain 

As  expected,  when  e  (which  is  directly  proportional  to  h)  approaches  zero 
this  equation  reduces  to  (10).  Equation  (19)  does  give  the  analytical  solution  to  the 
isospeed  case  but  is  overly  complicated  when  e  (thus  h)  is  small  (less  than  250  m). 
Equation  (10)  will  be  used  in  the  analysis  as  the  solution  for  Ar^z  and  G. 

D.  BILINEAR  WITH  MIXED  LAYER  SSP 

1,  Range  to  Outer  Edge  of  the  CZ 

With  the  addition  of  a  mixed  layer  to  the  bilinear  SSP,  and  remembering 
00  =  O',  the  initial  depressions  angle  8,  for  a  source  at  the  surface  is  expressed  by, 

02=acos 

2.  Range  to  the  Inner  Edge  of  the  CZ 

Figure  6  illustrates  the  SSP  with  the  addition  of  a  mbced  layer.  The  range 
equation  is  now 
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(20) 


Using  the  following  definitions 


H=1+!2 

h 


(21) 


n  = 


^2 

^0 


(22) 


(23) 


results  in  the  range  equation  becoming 


r=2R, 


rsinO,  -Hsin0n  +  — sin0, 
g2 


(24) 


r=2f? 

go 


'/l-Y'+e-  -H0o  +  !?/l-r72+e 

g2 


(25) 
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(25)  should  reduce  to  (9)  when  C2  is  equal  to  Cq,  and  gj  is  equal  to  gj. 
Using  these  conditions,  (25)  does  reduce  to  (9)  and  consequently  can  be  solved  by 
use  of  (10). 

Again,  finding  the  angle,  6q  =ein  at  which  r  =  rjnn„  , 


The  solution  for  Oj^  in  this  form  is  very  formidable.  However,  if  the 
surface  sound  speed  is  approximately  equal  to  the  sound  axis  speed,  then  «  /5“. 
Ojn  can  be  expressed  as 
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rj„n„  can  be  expressed  by: 


Which  reduces  to, 


Additionally,  rou,„  becomes: 


g2j 


(28) 
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3.  CZ  Gain 

The  width  of  the  CZ  is  now: 


Substitution  of  and  A  Fez  into  (2)  yields, 


Note  the  similarity  between  (13)  and  (30)  in  that  the  only  dependance  is 
upon  the  gradients. 

E.  COMPUTER  MODELS 

In  developing  the  computer  programs  to  validate  Oj^  and  G,  the  bilinear  SSP 
was  the  basic  case  of  interest  with  the  muced  layer  or  isospeed  layer  being  extensions. 
For  this  reason  several  conditions  are  placed  upon  the  computer  model  inputs. 

•  C2  is  always  less  than  Cq. 
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•  Cj  is  always  less  than  Cq. 

•  h  is  never  negative. 

•  The  mixed  layer  is  at  least  1  m. 

•  The  bottom  gradient  is  0.017  m/s/m. 

Appendices  B,  C,  and  D  contain  the  Fortran  listings  for  the  computer  programs 
used  to  validate  the  equations  for  Oj^  and  G  for  the  bilinear,  bilinear  with  isospeed 
layer,  and  bilinear  with  mixed  layer  SSFs  respectively.  Combinations  of  sound  speed, 
gradients,  and  thickness  of  isospeed  layer  were  used  where  appropriate.  But  in  order 
to  evaluate  the  different  cases,  a  range  for  the  input  values  needed  to  be  determined. 
Theses  ranges  were  obtained  from  historical  SSP’s  produced  by  ICAPS  and  contained 
in  [Ref.l],  Profiles  are  from  the  Pacific  (40N  MOW),  Atlantic  (3 IN  69W),  and 
Mediterranean  (36N  18E)  oceans  for  the  months  of  February,  May,  August,  and 
November  [Ref.  l:pp.l6-30].  Each  estimate  was  evaluated  using  the  following  ranges 
of  input  where  applicable: 

•  Co  ranges  from  1489  to  1541 

•  Cj  ranges  from  1475  to  1525 

•  Cj  ranges  from  1485  to  1525 

•  go  ranges  from  0.01  to  0.246 

•  gi  ranges  from  0.003  to  0.017 

•  g,  ranges  from  0.01  to  0.06 

•  h  ranges  from  0  to  250 
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These  sound  speed,  gradient,  and  h  ranges  were  programmed  into  the  validation 
programs,  and  equations  (10)/(13)  or  (26)/(30),  where  appropriate,  were  evaluated 
against  the  solution  ot),tained  by  ray  tracing  to  determine  Bjn ,  rjnn^,  Arc2  ,  and  G. 
The  input  data  range  increment  size  varied  depending  on  available  storage  space  for 
the  output  files.  As  the  number  of  variables  increased  the  increment  size  was 
adjusted  to  give  the  smallest  possible  value  while  maintaining  adequate  space  for  file 
.storage.  A  percent  error  for  G  was  calculated  using  the  following  algorithm: 

%  error  =  innv«tinrate-  ray  tracing 

ray  tracing 

The  historical  data  were  input  into  ICAPS  in  calculating  an  average  CZ  gain. 
[Ref.l:pp.75]  Appendix  E  contains  the  Fortran  listing  for  the  computer  program, 
CZGAIN,  used  to  compare  the  historical  data  from  ICAPS  to  the  results  for  G. 

Figure  7  depicts  a  typical  display  for  CZGAIN.  The  user  inputs  C2,  Cq,  do,  Ci,  dj, 
h,  and  bottom  depth.  Once  the  user  inputs  the  initial  values,  he  may  then 
manipulate  those  values  with  keyboard  entries.  Appends  F  contains  a  user’s  guide 
for  operation  of  CZGAIN.  All  of  the  data  are  color  coded  to  the  respective  SSP 
(which  labeled  1  and  2  for  the  black  and  white  depiction).  The  unit  dimensions  for 
the  Range  vs.  Theta  are  defined  at  the  top,  and  the  graph  is  scaled  so  that  each 
curve  starts  at  Theta(O),  Range(O).  As  Theta  (the  depression  angle  at  the  source) 
b  swept  from  Theta(O)  through  Theta(lO)  (0°  to  15°  as  shown  in  Fig.  7)  the  range 
at  which  that  ray  resurfaces  is  plotted.  Observe  in  Fig.  7  that  the  unit  dimensions 
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remain  the  same  for  both  curves,  and  each  curve  starts  from  the  same  point  on  the 
graph.  The  starting  reference  points  (Range(O)  and  Theta(O) )  are  also  color  coded 
to  their  respective  SSP’s.  Finally,  the  G’s  (also  color  coded)  are  presented.  The 
bilinear  approximation  for  rin„„  and  are  given  by  Range(+)  and  Theta(+). 
Notice  that  as  the  SSP’s  (1)  and  (2)  become  more  similar  the  Range  vs.  Theta  curves 
approach  each  other,  and  the  minimal  effect  from  the  isospeed  thickness  upon  the 
Range  vs.  Theta  curves. 

The  programs  contain  in  Appendix  B  and  C  contain  the  basic  algorithms  and 
equations  used  for  development  of  the  analytical  part  of  CZGAIN.  The  coding  for 
the  equations  may  more  easily  be  extracted  from  any  of  these  programs. 

All  of  the  programs  w«re  written  using  Microsoft’s  Fortran  Version  5.0  on  a 
Compaq  Deskpro  386/33  Model  84  computer.  The  programs  contained  in 
Appendices  B,  C,  and  D  have  also  run  on  the  Naval  Postgraduate  School’s 
mainframe  after  a  few  small  modifications  were  made  (e.g.  the  output  files  names 
were  modified).  A  special  software  patch  was  obtained  from  Microsoft  to  allow  the 
user  to  manipulate  the  data  parameters  without  having  to  use  the  "Enter"  key. 


Figure  7  -  CZGAIN  Display  Depiction 


IV.  RESULTS  AND  VALIDATION 


A.  GENERAL 

While  specific  results  will  be  presented  later,  overall  the  estimates  for  the  CZ 
gain  closely  approximated  the  values  obtained  by  ray  tracing  and  from  ICAPS. 

If  a  difference  of  more  than  3  dB  existed  between  the  estimate  and  the  solution 
from  ray  tracing  or  ICAPS,  the  estimate  was  considered  invalid.  The  3  dB  criterion 
was  chosen  because  the  numerous  estimations  to  real  ocean  conditions  made  not  only 
by  the  models  developed  here  but  also  by  ray  tracing.  A  tighter  requirement  upon 
the  models  serves  no  purpose  other  than  to  require  a  higher  degree  of  accuracy  from 
the  estimate,  forcing  additional  complexity  on  the  estimate. 

The  greatest  potential  for  range  errors  (and  thus  errors  in  the  gain  calculations) 
occurs  when  the  gradients  are  very  small  [Ref.l:p.66].  When  the  gradients  are  small 
an  extremely  high  degree  of  accuracy  is  required  for  6o  to  determine  .  The 
accuracy  is  required  to  diminish  the  effects  of  roundoff  error  and  decrease  the 
incrementation  size  for  Gq  .  As  Cj  and/or  Cj  approaches  Cq  the  gradients  approach 
zero  (assuming  that  the  respective  sound  speed  depths  remain  constant).  When  Cq 
and  Cl  and/or  Cq  and  c,  were  less  than  about  2  m/s  apart  the  model  becomes  invalid. 
Modifying  the  software  to  obtain  the  needed  degree  of  accuracy  for  Gq  was  not 
considered  to  be  necessary  as  this  closeness  of  sound  speeds  is  not  a  common 
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occurrence.  During,  validation,  the  sound  speeds  were  not  permitted  to  be  less  than 
2  m/s  apart. 

While  the  closeness  of  the  sound  speeds  rarely  occurs  in  the  ocean,  it  does 
demonstrate  a  weakness  in  the  model  which  would  not  have  been  discovered  without 
model  validation. 

In  comparing  this  model  to  the  results  from  earlier  models  [Ref.l:p.75]  it 
should  be  noted  that  earlier  results  used  a  source  and  receiver  at  18.3  m.  Also 
ICAPS  considers  reflection  from  the  surface,  absorption,  bottom  interactions,  and 
coherent  ray  interaction. 

The  tables  to  be  presented  contain  some  representative  samples  of  results  for 

•  G  from  (13)  or  (30) 

•  Ray  Tracing  (RT) 

•  Previous  calculator  models  (CALC)  [Ref.l:p.75] 

»  Averaged  ICAPS  values  [Ref.l:p.75j 

Not  all  of  the  data  are  presented  (the  entire  data  set  from  the  validation 
exceeds  2000  pages).  The  data  presented  are  a  subset  from  the  complete  date  set 
which  most  closely  matches  the  profiles  used  by  earlier  models.  The  tables  are 
intended  to  show  results  from  earlier  estimations  (CALC)  as  well  as  results  from 
more  complete  models  (ICAPS). 
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An  exact  fit  the  SSP’s  used  by  CALC  and  ICAPS  could  not  be  made  as  they  are 
linearized  by  different  methods.  In  determining  G  from  equations  (13)  and  (30)  the 
best  possible  fit  to  the  SSP’s  was  made. 

The  SSP’s  from  [Ref.l]  were  not  included  because  only  general  types  of  SSP 
and  not  specific  cases  of  profiles  were  not  examined.  The  SSP’s  used  are  readily 
available  in  [Ref.l]. 

The  ICAPS  data  is  an  averaged  value  and  is  obtained  from  [Ref.l]  by  the 
following  process: 

Two  levels  were  chosen  ...  which  bracket  the  majority  of  the  TL  points  within 
the  annulus.  The  approximate  midpoint  between  those  levels  was  then  picked 
as  "the"  TL  for  that  CZ.  [Ref.l:p.40] 

B.  BILINEAR  SSP 

Table  1  presents  the  data  resulting  from  use  of  equation  (13)  along  with  the 
RT,  CALC,  and  ICAPS  data.  The  SSP’s  numbered  18  through  20  [Ref. l:pp.  18-20] 
correspond  to  SSP’s  from  the  Pacific  Ocean  for  May,  August,  and  November 
respectively.  Note  the  agreement  between  the  results  from  equation  (13)  and  RT. 
Also  note  that  overall  a  closer  estimate  to  the  averaged  ICAPS  data  is  obtained  with 
equation  (13)  than  was  obtained  with  CALC. 
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TABLE  I  -  BILINEAR  SSP  CZ  GAINS  IN  DB  (PACIFIC  OCEAN) 


(13) 

RT 

ICAPS 

CALC 

18 

17.2 

17.2 

16.4 

15.8 

19 

18.2 

18.2 

15.8 

18.0 

20 

15.3 

15.3 

14.0 

10.0 

Table  II  presents  the  data  from  use  of  equation  (13)  along  with  the  RT,  CALC, 
and  ICAPS  data.  For  the  SSP’s  numbered  26  through  28  [Ref.l:pp.26-28] 
corresponding  to  the  SSP’s  from  the  Mediterranean  Sea  for  May,  August,  and 
November  respectively.  The  February  SSP  from  the  Mediterranean  Sea  did  not  have 
the  proper  gradient  structure  to  allow  the  formation  of  a  convergence  zone. 


TABLE  II  -  BILINEAR  SSP  CZ  GAINS  IN  DB  (MEDITERRANEAN  SEA) 


SSP  Number 

(13) 

RT 

ICAPS 

CALC 

26 

17.5 

17.7 

14.0 

18.0 

27 

19.6 

19.7 

17.7 

17.0 

28 

15.4 

15.4 

13.0 

17.0 

As  with  the  Pacific  Ocean  data,  the  results  from  equation  (13)  agree  closely 
with  ray  tracing  (less  than  1%  on  average).  Table  III  presents  the  percent  error 
between  the  averaged  ICAPS  data  and  equation  (13),  ray  tracing,  and  CALC. 
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TABLE  III  -  BILINEAR  SSP  -  AVERAGE  PERCENT  ERRORS  FROM  ICAPS 


(13) 

RT 

CALC 

13.9 

14.3 

18,2 

C.  BILINEAR  WITH  ISOSPEED  LAYER  SSP 

The  approximation  of  the  bilinear  case  (13)  was  used  rather  than  (19)  because 
of  the  complexity  involved  in  evaluating  0jn  by  (19).  Although  (19)  could  be 
programmed  on  a  computer  as  (13)  and  (30)  were,  the  thickness  of  the  isospeed  layer 
does  not  enter  significantly  into:the  range  equations  (unless  h  is  greater  than  250  m.). 
The  appropriate  SSP’s  [Ref.l:pp.21-24]  are  labeled  21  through  24  in  Table  IV.  These 
SSP’s  correspond  to  profiles  from  the  Atlantic  Ocean  for  February,  May,  August,  and 
November, 


TABLE  IV  .  BILINEAR  WITH  ISOSPEED  LAYER  SSP’S 
(ATLANTIC  OCEAN) 


SSP  Number 

(13) 

RT 

ICAPS 

CALC 

21 

16.7 

16.7 

15.5 

5.0 

22 

15.8 

15.7 

16.5 

16.0 

23 

17.4 

17.4 

16.0 

18.0 

24 

15.2 

15.2 

12.0 

17.0 

Again,  the  solution  obtained  from  use  of  (13)  agrees  closely  with  the  solution 
obtained  from  ray  tracing.  The  Atlantic  Ocean  SSP’s  were  more  difficult  to  segment 


because  of  a  more  complex  sound  speed  structure.  This  difficulty  has  also  been 
noted  in  previous  model  development  [Ref.l:p.66],  Table  V  presents  the  percent 
error  between  the  averaged  ICAPS  data  and  equation  (13),  ray  tracing,  and  CALC. 


TABLE  V  -  BILINEAR  WITH  ISOSPEED  LAYER  SSP  -  AVERAGE 
PERCENT  ERRORS  FROM  ICAPS 


(13) 

RT 

CALC 

11.9 

11.4 

31.2 

D.  BIUNEAR  WITH  A  MDCED  LAYER  SSP 

Table  VI  presents  the  data  resulting  from  use  of  equation  (30)  along  with  the 
ray  tracing,  ICAPS,  and  CALC  results.  Only  one  SSP,  numbered  17  [Ref.l:p.l7],  was 
appropriate  for  the  bilinear  with  mixed  layer  SSP  case.  This  SSP  corresponds  to  the 
SSP  from  the  Pacific  Ocean  for  February  and  was  the  only  profile  which  contained 
a  mixed  layer. 


TABLE  VI  -  BILINEAR  WITH  MIXED  LAYER  SSP’S 
(PACIFIC  OCEAN) 


SSP  Number 

(30) 

RT 

ICAPS 

CALC 

17 

17.3 

14.2 

13.4 

8.0 

Table  VII  present  the  percent  error  from  the  averaged  ICAPS  data  and 
equation  (30),  ray  tracing  and  CALC. 
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TABLE  VII  -  BILINEAR  WITH  MIXED  LAYER  SSP  - 
PERCENT  ERRORS  FROM  ICAPS 


wm 

RT 

CALC 

29.1 

6.1 

40.3 

The  unusually  high  percentage  error  for  (30)  in  this  case  is  a  direct  repercussion 
of  setting  the  surface  sound  speed  equal  to  the  sound  axis  speed.  In  fact  there  was 
about  a  12  m/s  difference  existed  between  C2  and  cl. 
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V.  CONCLUSIONS 


A.  MODEL  UMITATIONS 

Having  an  estimate  for  G  based  upon  a  bilinear  or  bilinear  with  mixed  layer 
SSP  would  be  sufficient  if  only  those  specific  SSP’s  existed.  Previous  models  have 
had  the  same  limitations  of  how  to  segment  the  SSP  [Ref.l:p.78].  A  balance  must 
be  maintained  between  practicality  and  desired  level  of  accuracy.  For  instance  in  this 
research  when  an  isospeed  layer  was  inserted  into  the  bilinear  SSP,  analysis  for 
quickly  became  unmanageable.  Notice  in  (30)  how  the  addition  of  a  simple  mixed 
layer  to  the  bilinear  SSP  added  to  the  complexity  of  (13).  The  addition  of  another 
layer  requires  two  more  parameters  to  be  introduced  and  will  increase  the  complexity 
of  the  gain  equation.  This  becomes  quite  complicated.  Though  higher  speed 
computers  (to  get  even  faster  in  the  near  future)  can  perform  complicated  algorithms 
swiftly  and  with  those  algorithms  can  come  even  closer  to  an  exact  solution  to  ray 
tracing,  the  intention  was  not  to  develop  an  exact  solution  but  a  good  estimate. 

This  model,  as  in  previous  work  [Ref.l:p.78],  does  not  check  for  the  bottom 
limited  sound  ray.  If  the  bottom  topography  interfered  with  part  of  the  beam  the 
Arcz  would  change  and  thus  G. 

As  noted  earlier,  with  the  sound  speeds  relatively  close  (less  than  about  2  m/s), 
the  model  becomes  invalid. 
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The  estimates  were  developed  for  a  source  and  receiver  on  the  surface. 
Consequently,  the  model  has  no  dependency  upon  source  or  receiver  depth. 

The  bottom  gradient  is  fixed  at  0.017  m/s/m.  Little  effort  would  be  required 
to  make  this  change,  but  the  operator  would  be  required  to  input  more  data. 
Additionally,  the  sound  speed  at  the  bottom  depth  may  not  be  known. 

CZGAIN  will  not  run  on  all  MS-DOS  machines.  High  resolution  graphics 
capability  is  required,  but  the  base  equations  (10),  (13),  (26),  and  (30)  can  be  coded 
into  practically  any  computer,  even  a  calculator  or  a  Macintosh. 

The  surface  sound  speed  and  sound  speed  at  the  sound  axis  must  be  equal  for 
(30)  to  be  valid.  Minor  deviations  from  this  requirement  may  still  be  evaluated  with 
(30),  but  the  extent  of  the  deviation  allowable  is  unknown. 

B.  MODEL  PRACTICAUTY 

In  determining  practicality,  the  desirable  qualities  presented  in  the  introduction 
will  be  considered. 

1.  Data  Requirements 

The  user/operator  is  required  to  have  an  SSP  simple  enough  so  that  no 
more  than  three  sound  speeds  and  three  depths  are  needed.  If  the  SSP  can  be 
approximated  by  a  bilinear  profile,  the  mbced  layer  information  may  be  omitted. 
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2.  Output 

The  output  of  CZGAIN  gives  the  user  a  repeat  of  the  inputs  in  text  and 
graphical  formats.  Additionally,  graphical  information  on  the  width  of  the  CZ  and 
range  to  the  inner  edge  of  the  CZ  are  presented. 

3.  Speed 

CZGAIN  runs  in  less  time  than  it  takes  to  input  the  data.  Any  of  the 
input  data  may  be  quickly  modified  without  having  to  reenter  all  data.  Earlier 
models  took  in  excess  of  nine  minutes  once  the  data  has  been  entered  [Ref.l:p.80]. 

4.  Program  Foundation 

Ray  tracing  theory  and  not  general  guidelines  or  thumbrules  were  used  in 
the  development  of  the  models. 

C.  COMPARISON  TO  EARLIER  MODELS 

The  results  from  the  model  agreed  closely  with  the  values  obtained  from  larger 
more  complicated  models  (ICAPS).  Furthermore,  the  estimates  for  G  developed 
were  closer  in  value  to  ICAPS  than  earlier  models  [Ref.l:p.75].  This  is  particularly 
interesting  since  earlier  models  used  a  more  highly  segmented  SSP  and  used  an 
iterative  solution  rather  than  an  analytical  one  [Ref.l:p.50].  Also,  a  tighter 
requirement  was  placed  on  this  model  for  an  acceptable  difference,  3  dB  vice  13  dB 
[Ref.l:p.81].  Validation  was  not  performed  on  the  programs  developed  previously 
[Ref.l]  so  no  validation  comparison  was  able  to  be  made. 
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D.  GENERAL 

The  estimates  for  CZ  gain  were  developed  using  ray  tracing  theory.  Equation 
(13)  is.  the  estimated  solution  for  G  for  the  bilinear  and  bilinear  with  isospeed  layer 
SSP’s.  Using  various  combinations  of  sound  speeds  and  gradients,  and  thickness  of 
the  isospeed  layer  where  appropriate,  the  solution  for  G  obtained  from  (13)  agreed 
more  closely  with  the  solution  for  G  from  ray  tracing  as  well  as  ICAPS.  Equation 
(30)  is  the  estimated  solution  for  G  for  the  bilinear  with  mixed  layer  SSP.  This 
equation  was  developed  under  the  condition  of  the  surface  and  sound  axis  sound 
speeds  being  equal.  As  long  as  this  condition  was  met,  the  solution  from  (30)  and 
solution  from  ray  tracing  agreed  fairly  well  (no  ICAPS  data  was  available).  As  long 
as  the  SSP’s  do  vary  from  this  requirement,  equations  (13)  and  (30)  give  a  reasonable 
estimates  of  the  convergence  zone  gain. 
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APPENDIX  A  r  SERIES  EXPANSIONS 


Several  expansions  are  used  throughout  the  development  of  various  equations. 
The  following  is  a  synopsis  of  those  used. 


A.  TRIGONOMETRIC  SERIES  EXPANSIONS 
1.  Cosine  0 


cos  0  =  1-— +  —  ■ 

2!  4!  6!  8! 


2.  Tangent  0 


t.nfl  flx0^20\l70’  620'" 

tan  0=0  +  —  +  —  + - + - + 

3  15  315  2835 


where  |0| 


B.  BINOMIAL  SERIES  EXPANSION 


\2 


1-0 


=  1+20+30- +403+50^+605  +  « 


where  0”  <  1 
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APPENDIX  B 


CC  RESWPT2.F0R 
CC 

CG  Program  to  evaluate  the  estimation  angle  to  the  inner  edge  of 
CC  the  reswept  zone  and  the  CZ  gain.  An  iteration  of  the  range 
CC  equation  is  performed  to  obtained  the  true  value.  Finally,  the 
CC  percent  error  of  the  angle  of  the  ray  to  the  inner  edge  and  CZ  gain 
CC  as  compared  to  the  iteration  value  are  calculated. 

CC 

CC  by  A1  Jenkins 
CC 

REAL*8  c_0,  c_l,  g_0,  g_l,  R_0,  R_l,  theta_0,  theta_l, 

+  range_e_min,  rangej_min,  temp_range,  G,  gamma, 

+  range_outer,  theta_e,  thetaj, 

+  Gainjter,  Gain_est 

LOGICAL  theta__mi.n 
INTEGER  i,  counter 

OPEN  (UNIT=1,  ACCESS  =  ’APPEND’,  FILE=  ’reswpt2a.dat’, 

+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=2,  ACCESS^  ’APPEND’,  FILE=  ’reswpt2b.dat’, 

+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=3,  ACCESS  =  ’APPEND’,  FILE=  ’reswpt2c.dat’, 

+  STATUS  =  ’UNKNOWN’) 

WRITE(1,9000) 

theta_min  =  .FALSE, 
counter  =  0 

g  1  =  0.003 
DO  400  m  =  1,  10 
WRITE(  *,•  )  m 
c_l  =  1475 
DO  300  1  =  1,  10 
g  0  =  .017 
DO  200  k  =  1,  10 
c  0  =  1489 
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DO  100  j  =  1,10 

IF(  c_0  .GT.  c_l )  THEN 


G  =  1  +  g_0  /  g_l 
gamma  =  c_l  /  c_0 

theta_min  =  .FALSE. 

theta_0  =  SQRT(  (l-gamma**2)  /  (G**2-l) ) 
theta_l  =  acos(  c_l  •  cos(theta_0)  /  c_0  ) 
theta_e  =  theta_0 

R_0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R”1  =  c_l  /  (  g_l  *  cos(theta_l) ) 

range_e_min  =  R_0  *  (  sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l) 

theta_0  =  0 

theta_l  =  acos(  c_l  *  cos(theta_0)  /  c_0  ) 

=  c__0  /  (  g_0  *  cos(theta_0) ) 

R~1  =  c_l  /  (  g_l  •  cos(theta_l) ) 

range_outer  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l) 

temp__range  =  range_outer 

DO  10  i  =  1,10000 
IF(  .NOT.  theta_min  )  THEN 
theta_0  =  theta_0  -r  0.000035 
theta_l  =  acos(  c_l  •  cos(theta_0)  /  c_0  ) 

R_0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l  •  cos(theta_l) ) 

range_i_min  =  R_0  *  (  sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l) 

CC  Loops  set  up  so  as  to  write  only  the  minimum  angle  and  range 
IF  (range_i_min.LT.temp_range)  THEN 
temp_range = range_i_min 
ELSE  IF  (.NOT.  theta_min)  THEN 
range_i_min  =  temp_range 
theta_min = .TRUE, 
thetaj  =  theta  0 
ENDIF 
ENDIF 
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10  CONTINUE 


error  =  (theta_i  -  theta_e)  /  thetaj 
error  =  100  *  ^S(  error) 

WRITE(2,9010)  error 
IF(  error  .GE.  5  ) 

+  WRITE(1,9001)  c_0,  g_0,  c_l,  g_l,  theta_e,  thetaj,  error 

counter  =  counter  +  1 

Gainjter  =  10  *  LOG10(  2*  range J_min  / 

+  (  range  outer  -  range  i  min) ) 

Gain^est  =  fO  *  LOG10(  2*“SQRT(  G**2  -  1)  /  (G  -  SQRT(G**2  -1))) 
error  =  (Gainjter  -  Gain_est)  /  Gainjter 
error  =  100  *  ABS(  error) 

WRITE(3,9010)  error 

IF(  (error*GainJter)  .GE.  3  ) 

+  WRITE(1,9002)  Gain  iter,  Gain  est,  theta_e,  theta  i,  error 

END  IF 

c  0  =  c_0  +  6 
100  CONTINUE 

g  0  =  g  0  +  0.0229 
200  CONTINUE 
c  1  =  c  1  +  5 
300  contTnue 

g  1  =  g  1  +  0.002 

400  CONTINUE 

WRITE(  1,  •  )  counter 
CLOSE  (1,  STATUS  =  ’KEEP) 

CLOSE  (2,  STATUS  =  ’KEEP’) 

9000  FORMAT  (lX,’c0’,8X,’g0’,8X,’cl’,8X,’gl’,8X,’t-e’,6X,’t-i’,6X, 

+  ’%  Error’) 

9001  FORMAT  (1X,F6.1,3X,F5.3,6X,F5.0,4X,F5.3,6X,F5.4,4X,F5.4,4X,F9.6) 

9002  FORMAT  (1X,F6.2,3X,F6.2,6X,F5.4,4X,F5.4,4X,F9.6) 

9010  FORMAT  (1X,F13.8) 

9100  FORMAT(  IX,  F7.3,  4X,  F7.3,  3X,  F6.1,  3X,  F5.3,  6X,  F5.0,  4X, 

+  F5.3,  2X,  F4.1  ) 

9101  FORMAT(  IX,  F4.1,  4X,  F4.1) 

END 
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APPENDIX  C 

CC  RESWPT3.F0R 
CC 

CC  Program  to  evaluate  the  estimation  angle  to  the  inner  edge  of 
CC  the  reswept  zone  and  the  CZ  gain.  An  iteration  of  the  range 
CC  equation  is  performed  to  obtained  the  true  value.  Finally,  the 
CC  percent  error  of  the  angle  of  the  ray  to  the  inner  edge  and  CZ  gain 
CC  as  compared  to  the  iteration  value  are  calculated,  h  is  the  thickness 
CC  of  the  isospeed  layer 
CC 

CC  by  A1  Jenkins 
CC 

REAL’S  c_0,  c_l,  g_0,  g_l,  R_0,  R_l,  theta_0,  theta_l, 

+  range_e_min,  rangej_min,  temp__range,  G,  gamma, 

+  range__outer,  theta_e,lhetaj, 

+  Gainjter,  Gain_est,  h 

LOGICAL  theta_min 
INTEGER  i,  j,  k,  1,  m,  n,  counter 

OPEN  (UNIT=1,  ACCESS=  ’APPEND’,  FILE=  ’reswpt3a.dat’, 

+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=2,  ACCESS^  ’APPEND’,  FILE=  ’reswpt3b.dat’, 

+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=3,  ACCESS=  ’APPEND’,  FILE=  ’reswpt3c.dat’, 

+  STATUS  =  ’UNKNOWN’) 

WRITE(1,9000) 

theta_min  =  .FALSE, 
counter  =  0 
h=0 

DO  500  n  =  1,10 
g_l  =  0.003 
DO  400  m  =  1,  10 
WRITE(  )  n,m 
c  1  =  1475 
DO  300  1  =  1,  10 
g_0  =  .017 
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DO  200  k  =  1,  10 
c  0  =  1489 
DO  100  j  =  140 

IF(  c_0  .GT.  c_l )  THEN 

G  =  1  +  g_0  /  g_l 
gamma  =  c_l  /  c_0 

theta_min  =  .FALSE. 

theta_0  =  SQRT(  (l-gamma**2)  /  (G**2-l) ) 
theta_l  =  acos(  c_l  •  cos(theta_0)  /  c_0  ) 
theta_e  =  theta_0 

R_0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l  *  cos(theta_l) ) 

range_^e_min  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  h  *  cotan(  theta_l ) 

theta_0  =  0 

theta_l  =  acos(  c_l  *  cos(theta_0)  /  c_0  ) 

R_0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l  •  cos(theta_l) ) 

range_outer  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  h  *  cotan(  theta_l  ) 

temp_range  =  range_outer 

DO  10  i  =  1,10000 
IF(  .NOT.  theta_min  )  THEN 
theta_0  =  theta_0  +  0.000035 
theta_l  =  acos(  c_l  *  cos(theta_0)  /  c_0  ) 

R__0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l  *  cos(theta_l) ) 

range_i_min  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  h  *  cotan(  theta_l  ) 

CC  Loops  set  up  so  as  to  write  only  the  minimum  angle  and  range 
IF  (range_i_min.LT.temp_range)  THEN 
temp  range =range_i_min 
ELSETF  (.not.  theta_min)  THEN 
range_i_min  =  temp_range 
theta_min=.TRUE. 
theta  i  =  theta  0 
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ENDIF 

ENDIF 

10  CONTINUE 

error  =  (theta_i  -  theta_e)  /  theta_i 
error  =  100  *  ^S(  error) 

WRITE(2,9010)  error 
IF(  error  .GE.  5  ) 

+  WRITE(1,9001)  c_0,  g_0,  c_l,  g_l,  theta_e,  theta_i,  error 

counter  =  counter  +  1 

Gainjter  =  10  *  LOG10(  2*  range_i_min  / 

+  (  range  outer  -  range  i  min)  3 

Gain_est  =  fO  *  LOG10(  2*“SQRT(  G**2  •  1)  /  (G  -  SQRT(G**2  -1))) 
error  =  (Gainjter  -  Gain_est)  /  Gainjter 
error  =  100  •  ABS(  error) 

WRITE(3,9010)  error 
IF(  (error’Gain  iter)..GE.  3  )  THEN 
WRITE(l,900l)  c_0,  g_0,  c_l,  g_l,  1. 

WRITE(1,9002)  Gain  iter,  Gain  est,  theta  e,  theta j,  error 
ENDIF 
END  IF 

c  0  =  c  0  +  5.2 
100  CONTINUE 
g  0  =  g  0  +  0.0229 
200  CONTINUE 
c  1  =  c  1  +  5 
300  CONTINUE 

g_l  =  g_l  +  0.002 
400  CONTINUE 
h  =  h  +  25 
500  CONTINUE 

WRITE(  1,  •  )  counter 
CLOSE  (1,  STATUS  =  ’KEEP’) 

CLOSE  (2,  STATUS  =  ’KEEP’) 

9000  FORMAT  (IX, ’cO’,8X,’gO’,8X, ’cl’, 8X,’gl’,8X,’t-e’,6X,’t-i’.6X, 

+  ’%  Error’) 

9001  FORMAT  (IX, F6.1,3X,F5.3,6X,F5.0,4X,F5.3,6X,F5.1,4X,F5.4,4X,F9.6) 

9002  FORMAT  (1X,F6.2,3X,F6.2,6X,F5.4,4X,F5.4,4X,F9.6) 

9010  FORMAT  (1X,F13.8) 

9100  FORMAT(  IX,  F7.3,  4X,  F7.3,  3X,  F6.1,  3X,  F5.3,  6X,  F5.0,  4X, 
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+  F5.3,  2X,  F4.1 ) 

9101  FORMAT(  IX,  F4.1,  4X,  F4.1) 
END 
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APPENDIX  D 


CC  RESWPT4.F0R 
CC 

CC  Program  to  test  the  estimated  value  for  theta  against  the 
CC  iterative  value  of  theta  for  the  inner  edge  of  the  reswept  zone 
CC  of  a  convergence  zone. 

CC 

CC  by  A1  Jenkins 
CC 


REAL*8  c_0,  c_l,  g_0,  g_l,  R_0,  R_l,  theta_0,  theta_l, 

+  range_e_min,  range_i_min,  temp_range,  G,  gamma, 

+  range_outer,  theta_e,  thetaj,  height,  c_2,  g_2, 

+  H,  eta,  R_2 

LOGICAL  theta_min 
INTEGER  i,  j,  k,  1,  m,  n,  o,  p,  counter 

OPEN  (UNIT=1,  ACCESS  =  ’APPEND’,  FILE=  ’RESWPT4A.DAr, 
+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=2,  ACCESS=  ’APPEND’,  FILE=  ’RESWPT4B.DAr, 
+  STATUS  =  ’UNKNOWN’) 

OPEN  (UNIT=3,  ACCESS  =  ’APPEND’,  FILE=  ’RESWPT4CDAr, 
+  STATUS  =  ’UNKNOWN’) 

WRITE(1,9000) 

theta_min  =  .FALSE, 
counter  =  0 


g_2  =  0.01 

DO  700  p  =  1,10 

c  2  =  1475 
DO  600  0  =  1,10 

height  =  0 
DO  500  n  =  1,1 
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g_l  =  0.008 
DO  400  m  i=  1,  5 
WRITE(  *,*  )  TUi  n,  0,  p 


c  1  =  1475 
DO  300  1  =  1,  1 


g  0  =  .01 

DO  200  k  =  1,  1 

c  0  =  1476 
DO  100  j  =  1,10 

IF(  c_0  .GT.  c_l  .AND.  c_0  .GT.  c_2)  THEN 

G  =  1  +  g_0  /  g_l 
gamma  =  c_l  /  c_0 
H  =  1  +  g_0  /  g_2 
eta  =  c_2  /  c_0 

theta_min  =  .FALSE. 

theta  0  =  H  *  SQRT(  1  -  gamma* *2)  /  SQRT(  (  G  +  g  0/g  2)* *2  - 
+  "H**2 ) 

theta_l  =  acos(  gamma  *  cos(theta__0) ) 
iheta_2  =  acos(  eta  *  cos(theta_0)  jT 
theta_e  =  theta_0 

R_0  =  c_0  /  (  g_0  *  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l:  *  cos(theta_l) ) 

R_2  =  c_2  /  (  g_;i  *  cos(theta_2) ) 

range_e_min  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  R_2  *  (  sin(theta_2)  -  sin(theta_0) ) 

theta_0  =  0 

theta_l  =  acos(  gamma  *  cos(theta_0) ) 
theta_2  =  acos(  eta  *  cos(theta_0) ) 

R_0  =  c_0  /  (  g_0  *  cos(the7a_0) ) 

R_1  =  c_l  /  (  g_l  *  cos(theta_l) ) 

R__2  =  c_2  /  (  g_2  *  cos(theta_2) ) 

range_outer  =  R_0  *  ( sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  R_2  *  (  sin(theta_2)  -  sin(theta_0) ) 

temp_range  =  range_outer 
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DO  10  i=  1,10000 
iF(  .NOT.  theta_min  )  THEN 
theta_0  =  theta_0  +  0.000035 
theta_l  =  acos(  gamma  *  cos(theta_0) ) 
theta_2  =  acos(  eta  *  cos(theta_0) ) 

R_0  =  c_0  /  (  g_0  •  cos(theta_0) ) 

R_1  =  c_l  /  (  g_l  *  cos(theta_l) ) 

R_2  =  c_2  /  (  g_2  •  cos(theta_2) ) 

range J_min  =  R_0  *  (  sin(theta_l)  -  sin(theta_0) )  + 

+  R_1  *  sin(theta_l)  +  R_2  *  ( sin(theta_2)  -  sin(theta_0) ) 

CC  Loops  seF  up  so  as  to  write  only  the  minimum  angle  and  range 
IF  (range_i_min.LTtemp_range)  THEN 
temp  range = range  i_min 
ELSE  IF  (.NOT.  thera_min)  THEN 
theta_min = .TRUE, 
thetaj  =  theta_0 
ENDIF 
ENDIF 

10  CONTINUE 

IF(  .NOT.  theta_min  )  WRITE(  1,*  )  counter 
error  =  (thetaj  -  theta_e)  /  thetaj 
error  =  100  *  ^S(  error) 

WRITE(2,9010)  error 
IF(  error  .GE.  5  ) 

+  WRITE(  1,9001)  c_0,  g_0,  c_l,  g_l,  theta_e,  height,  error 
counter  =  counter  +  1 

Gainjter  =  10*LOG10(2*range  i_min/(range_outer*rangej_min) 
Gain_est  =  10*LOG10(2*SQRT((G+g_0/g  2)**2-H**2)/((G+g_0/g_2)- 
+  SQRT((g+g_0/g_2)**2-H**2)) 

error  =  (Gainjter  -  Gain_est)  /  Gainjter 
error  =  100  *  ABS(  error) 

WRITE(3,9010)  error 
IF(  (error* Gainjter)  .GE.  3  )  THEN 
WRITE(  1,9001)  c_0,  g_0,  c_l,  g_l,  height 
WRITE(  1,9002)  Gainjter,  Gain_est,  theta_e,  thetaj,  error 
ENDIF 

ENDIF 


c_0  =  c_0  +  7.5 
100  CONTINUE 
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g_,0  =  g_0  +  0.04 
200  CONTINUE 

c_l  =  c_l  +  5 
300  CONTINUE 


g_l  =  g_l  +  0.004 
400  CONTINUE 

height  =  height 
500  CONTINUE 

c  2  =  c  2  +  5 
600  CONTINUE 

g  2  =  g_2  +  0.005 
700  CONTINUE 

WRITE(  1,  *  )  counter 
CLOSE  (1,  STATUS  =  ’KEEP’) 

CLOSE  (2,  STATUS  =  ’KEEP’) 

9000  FORMAT  (IX, ’cO’,8X,’gO’,8X, ’cl’, 8X,’gl’,8X,’t-e’,6X,’h’,6X, 

+  ’%  Error’) 

9001  FORMAT  (IX, F6.1,3X,F5.3,6X,F5.0,4X,F5.3,6X,F5.4,4X,F5.0,4X,F13.6) 

9002  FORMAT  (1X,F6.2,3X,F6.2,6X,F5.4,4X,F5.4,6X,F13.6) 

9010  FORMAT  (1X,F13.8) 

END 
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APPENDIX  E 


CC  GAIN.FOR 
CC 

CC  Program  for  finding  the  CZ  gain 
CC 

CC  by  A1  Jenkins 
CC 

INCLUDE  TGRAPH.FF 
INCLUDE  TGRAPH.FD’ 

LOGICAL  fourcolors 
EXTERNAL  fourcolors 
REALMS  cO,  cl,  c2,  h,  dO,  dl,  d3 

c  WRITE  (*,*)  ’  Enter  cO  -  ’ 
c  READ  (*,*)  cO 
c  WRITE  (*,*)  ’  Enter  dO  -  ’ 
c  READ  (*,*)  dO 
c  WRITE  (•,•)  ’  Enter  cl  -  ’ 
c  READ  (•,*)  cl 
c  WRITE  (*,*)  ’  Enter  dl  -  ’ 
c  READ  (*,•)  dl 
c  WRITE  (*,•)  ’  Enter  c2  -  ’ 
c  READ  (*,*)  c2 
c  WRITE  (*,•)  ’  Enter  h  -  ’ 
c  READ  (*,*)  h 
c  WRITE  (*,*)  ’  Enter  d3  -  ’ 
c  READ  (*,*)  d3 

cO  =  1500 
cl  =  1475 
c2  =  1490 
h  =  250 
dO  =  100 
dl  =  500 
d3  =  3000 
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IF  (  fourcolorsO  )  THEN 

CALL  plot  all(cO,  cl,  c2,  h,  dO,  dl,  d3) 

ELSE 

WRITE  (*,*)  ’  This  program  requires...’ 

END  IF 
END 
CC 

cc 

CC  FOURCOLORS- 
CC  p.180  Advanced  Topics  Guide 
CC  Microsoft  Corporation 
CC 

LOGICAL  FUNCTION  fourcolors() 

INCLUDE  ’FGRAPH.FD’ 

INTEGER*2  dummy 
RECORD  /videoconfig/  screen 
COMMON  screen 

C 

C  Set  to  max  number  of  available  colors. 

C 

CALL  getvideoconfig(  screen  ) 

SELECT  CASE(  screen.adapter ) 

CASE(  $CGA,  $OCGA  ) 
dummy  =  setvideomode(  $MRES4COLOR  ) 
CASE(  $EGA,  $OEGA  ) 
dummy  =  setvideomode(  SERESCOLOR  ) 
CASE(  SVGA,  SOVGA ) 
dummy  =  setvideomode(  $VRES16COLOR  ) 
CASE  DEFAULT 
dummy  =  0 
END  SELECT 

CALL  getvideoconfig(  screen  ) 
fourcolors  =  .TRUE. 

IF(  dummy  .EQ.  0  )  fourcolors  =  .FALSE. 

END 

CC 

CC 
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CC  PL0T_ALL  Displays  the  graphs 

CC 

CC 

SUBROUTINE  plot_all(c0,  cl,  c2,  h,  dO,  dl,  d3) 

INCLUDE  TGRAPH.FD’ 

REALMS  cO,  cl,  c2,  h,  dO,  dl,  d3 

INTEGER*2  dummy,  halfx,  halfy,  xwidth,  yheight 

INTEGER  *2  cols,  rows 

RECORD  /videoconfig/  screen 

COMMON  screen 

CHARACTER*!  string,  getch 

RECORD  /rccoord/  curpos 

10  CONTINUE 

C 

C  Find  screen  dimensions 
C 

CALL  clearscreen  (  SGCLEARSCREEN  ) 

xwidth  =  screen.numxpixels 

yheight  =  screen.numypixels 

cols  =  screen.numtextcols 

rows  =  screen.numtextrows 

halfx  ==  xwidth/2  *  0.8 

halfy  =  (yheight/rows)*  (rows/2) 

C 

C  Sound  Speed  window 
C 

CALL  setviewport(  0,  0,  halfx  ,  halfy  ) 

CALL  settextwindow(  1,  1,  rows  /  2,  cols  /  2  ) 
dummy  =  setwindow(  .TRUE.,  -2.0,  -1.2,  2.0, 1.75  ) 

CALL  sound_speed(c0,  cl,  c2,  h,  dO,  dl,  d3) 
dummy  =  rectangle  (SGBORDER,  0,  0,  halfx  -  4,  halfy  -  2  ) 
C 

C  Range  vs.  Theta  window 
C 

CALL  setviewport(  halfx  -  2,  0,  xwidth-1,  yheight  - 1  ) 

CALL  settextwindow(  1,  (cols/2-5) ,  rows,  cols ) 
dummy  =  setwindow(  .TRUE.,  -2.0,  -2.0,  2.0,  2.0  ) 

CALL  rangej)lot(c0,  cl,  c2,  h,  dO,  dl) 

dummy  =  rectangle_w(  SGBORDER,  -2.0,  -2.0,  2.0,  2.0  ) 
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ooo  non 


Graph  Description  Window 

CALL  setviewport(  0,  halfy+1 ,  halfx-1  ,  yheight-1) 

CALL  settextwindow(  rows  /  2  +  1,  0,  rows,  cols  /  2  ) 
dummy  =  setwindow(  .TRUE.,  -2.0,  -1.2,  2.0, 1.75  ) 

CALL  graph_description(  cO,  dO,  cl,  dl,  d3,  h,  c2  ) 
dummy  =  rectangle  ($GBORDER,  0,  0,  halfx  -  4,  halfy  -  2  ) 


Check  for  increasing,  adding,  rasing  sound  speed  at  surface 


CALL  settextposition(  0,  0,  curpos ) 
string  *  getchO 
SELECT  CASE(  string  ) 

CASE(  T’ ) 

GOTO  20 
CASE(  ’!’,  ’A’,  ’a’ ) 

GOTO  25 
CASE(  ’2’ ) 

GOTO  30 
CASE(  ’S’,  ’s’ ) 

GOTO  35 
CASE(  ’3’ ) 

GOTO  40 
CASE(  ’#’,  ’D’,  ’d’ ) 

GOTO  45 
CASE(  ’4’ ) 

GOTO  50 
CASE(  ’$’,  ’F,  ’f  ) 

GOTO  55 
CASE(  ’5’ ) 

GOTO  60 

CASE(  ’G’,  ’g’ ) 

GOTO  65 
CASE(  ’6’ ) 

GOTO  70 
CASE(  ’H’,  ’h’  ) 

GOTO  75 

CASE(  ’E’,  ’e’,  ’Q’,  ’q’) 

GOTO  100 
CASE  DEFAULT 
GOTO  10 
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END  SELECT 
GOTO  10 


20  c2  =  c2  -  1 

IF  (c2.LE.cl)  c2  =  c2  +  1 
IF  (c2.LE.O  )  c2  =  c2  +  1 
GOTO  10 
25  c2  =  c2  +  1 

IF  (c2.GE.cO)  c2  =  c2  -  1 
GOTO  10 
30  cO  =  cO  - 1 

IF  (cO.LE.c2)  cO  =  cO  +  1 
IF  (cO.LE.cl)  cO  =  cO  +  1 
IF  (cO.LE.O  )  cO  =  cO  +  1 
GOTO  10 
35  cO  =  cO  +  1 
GOTO  10 
40  dO  =  dO  .  5 

IF  (dO.LE.O)  dO  =  dO  +  5 
GOTO  10 
45  dO  =  dO  +  5 

IF  (dO.GE.dl)  dO  =  dO  -  5 
GOTO  10 
50  cl  =  cl  •  1 

IF  (cl.LE.O  )  cl  =  cl  +  1 
GOTO  10 
55  cl  =  cl  +  1 

IF  (cl.GE.c2)  cl  =  cl  - 1 
IF  (cl.GE.cO)  cl  =  cl  -  1 
GOTO  10 
60  dl  =  dl  -  5 

IF(  dl  .LE.  dO  )  dl  =  dl  +  5 
GOTO  10 
65  dl  =  dl  +  5 

IF  (dl.GE.(d3  -  h))  dl  =  dl  -  5 
GOTO  10 
70  h  =  h  -  5 

IF(  h  .LT.  0  )  h  =  h  +  5 
GOTO  10 
75  h  =  h  +  5 

IF(  h  .GE.  (d3  -  dl  -  dO) )  h  =  h  -  5 
GOTO  10 
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100  CONTINUE 

dummy  =  setvideomode  (  SDEFAULTMODE  ) 
END 


CC 

CC 

CC  SOUND  SPEED  Subroutine  to  plot  the  sound  speed  profile 

CC 

CC 

SUBROUTINE  sound jpeed(c0,  cl,  c2,  h,  dO,  dl,  d3) 

INCLUDE  TGRAPH.FD’ 

REAL’S  cO,  cl,  c2,  h,  dO,  dl,  d3 
INTEGER*2  dummy 
DOUBLE  PRECISION  x,  y,  scale_x 
RECORD  /videoconfig/  screen 
RECORD  /wxycoord/  wxy 
RECORD  /rccoord/  curpos 
COMMON  screen 
CHARACTER’S  string 

dummy  =  settextcolor(  15  ) 

CALL  settextposition(  8,  2,  curpos ) 

CALL  outtext(  ’Depth’ ) 

CALL  settextposition(  2,  10,  curpos ) 

CALL  outtext(  ’Sound  Speed’ ) 

CALL  settextposition(  3,  12,  curpos ) 
dummy  =  settextcolor(  14  ) 

WRITE  (string,  ’(I4)’)INT(  c2  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  3,  19,  curpos  ) 
dummy  =  settextcolor(  12  ) 

WRITE  (string,  ’(I4)’)INT(  cO  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

C 

C  Draw  border 
C 

dummy  =  setcolor(  1  ) 

dummy  =  rectangle_w(  SGBORDER,  -1.00,  -1.00,  1.00,  1.00  ) 
dummy  =  rectangle_w(  SGBORDER,  -1.02,  -1.02,  1.02,  1.02  ) 


61 


o  o  o  o  n  o 


Plot  points 

scalejc  =  -  50.0 
dummy  =  setcolor(  14  ) 

CALL  moveto_w(  0, 1,  wxy  ) 

X  =  scale  X  •  (  1  -  cO  /  c2  ) 
y  =  1  -  2"*  dO  /  d3 
dummy  =  lineto_w(  x,  y  ) 

CALL  moveto_w(  x,  y,  wxy  ) 

X  =  scale  X  *  ( 1  *  cl  /  c2  ) 
y  =  1  -  2“*  dl  /  d3 
dummy  =  lineto_w(  x,  y  ) 

CALL  moveto  w(  x,  y,  wxy  ) 
y  =  1  -  2  *  (  dl  +  h  )  /  d3 
dummy  =  lineto_w(  x,  y  ) 

CALL  moveto_w(  x,  y,  wxy  ) 

X  =  scale_x  •(l-(cl  +  0.017*(d3-dl-h))/c2) 
dummy  =  lineto_^w(  x,  -1 ) 
dummy  =  setcolor(  3  ) 

Plot  points  for  bilinear  profile 

dummy  =  setcolor(  12  ) 

X  =  scale^x  *  (  1  -  cO  /  c2  ) 

CALL  moveto_w(  x,  1,  wxy  ) 
x  =  scale  X  *  ( 1  -  cl  /  cO  ) 
y  =  1  -  2“*  dl  /  d3 
dummy  =  lineto_w(  x,  y  ) 

CALL  moveto  w(  x,  y,  wxy  ) 

X  =  scale_x  *(l-(cl+0.017*(d3-dl))/c0) 
dummy  =  lineto_w(  x,  -1 ) 
dummy  =  setcoIor(  3  ) 

END 
CC 
CC 

CC  RANGE  PLOT  Subroutine  for  plotting  range  versus  theta 

CC 

CC 

SUBROUTINE  range J3lot(c0,  cl,  c2,  h,  dO,  dl) 
INCLUDE  ’FGRAPH.FD’ 
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•INTEGER*2  dummy,  i 

DOUBLE  PRECISION  x,  y,  theta_2,  rangejcale 

DOUBLE  PRECISION  theta,  theta_graph,  range,  range_graph 

REAL*8  unit_dimension,  range_12,  theta_12,  gamma,  G 

REAL’S  range_14,  theta_14 

REAL’S  cO,  cl,  c2,  dO,  dl,  h,  scale_r,  scale_t 

RECORD  /videoconfig/  screen 

RECORD  /wxycoord/  wxy 

RECORD  /rccoord/  curpos 

COMMON  screen 

CHARACTER’S  string 

Label  grid 

dummy  =  settextcolor(  13  ) 

CALL  settextposition(  14,  0,  curpos ) 

CALL  outtext(  ’Range’ ) 

CALL  settextposition(  15,  8,  curpos ) 

CALL  outtext(  ’  0’ ) 

CALL  settextposition(  10,  8,  curpos ) 

CALL  outtext(  ’  2’ ) 

CALL  settextposition(  4,  8,  curpos ) 

CALL  outtext(  ’  4’ ) 

CALL  settextposition(  21,  8,  curpos ) 

CALL  outtext(  ’-2’ ) 

CALL  settextposition(  27,  8,  curpos ) 

CALL  outtext(  ’-4’ ) 

CALL  settextposition(  35,  11,  curpos ) 

CALL  outtext(  ’0  2  4  6  8  Theta’ ) 


Draw  border 
dummy  =  setcolor(  1  ) 

dummy  =  rectangle_w(  SGBORDER,  -1.00,  -1.80,  1.90,  1.90  ) 
dummy  =  rectangle_w(  SGBORDER,  -1.02,  -1.82,  1.92,  1.92  ) 

Draw  Grid 


X  =  -0.71 
y  =  -1.48 
DO  i  =  1,  9 
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dummy  =  setcolor(  8  ) 

CALL  moveto_w(  x,  -1.8,  wxy  ) 
dummy  =  lineto_w(  x,  1.9  ) 

CALL  moveto_w(  -1.0,  y,  wxy  ) 
dummy  =  lineto_w(  1.9,  y  ) 
y  -  y  +  0.37 

x  =  X  +  0.29 

END  DO 
C 

C  Plot  curve 
C 

gO  =  ABS(  (cO  -  cl)  /  (dl  -  dO) ) 

gl  =  ABS(  0.017  ) 

g2  =  ABS(  (cl  -  c2)  /  dO  ) 

C  Set  scaling  factor  for  vertical  axis,  based  on  bilinear  profile 
gamma  =  cl  /  cO 
G  =  1  +  gO  /  gl 

theta  =  SQRT(  (  1  -  gamma**2  )  /  (  G**2  - 1 ) ) 
theta_l  =  acos  (cl  *  cos(theta)  /  cO) 

theta_2  =  acos  (c2  *  cos(theta)  /  cO) 

c  range=(cO/(gO*cos(theta)))*((l+gO/gl)*DSQRT(l-(cl/cO)**2  * 
c  +  cos(theta)**2)  -  sin(theta))  + 

c  +  h*cotan(acos((cl/c())*cos(theta))) 

range  =  c2/(g2*cos(theta_2))  *  (sin(theta_2)  -  sin(theta))  + 

+  c0/(g0*cos(theta))  *  (sin(theta_l)  -  sin(theta))  + 

+  h  *  cotan(theta_l)  + 

+  cl/(gl*cos(theta_l))  *  sin(theta_l) 

range_scale  =  0.045  *  range 

dummy  =  setcolor(  14  ) 
scale_r  =1.0 

C  Set  horizontal  scale  so  that  1.5  degrees  equals  0.1  across  the 
C  total  width  of  2.9 

scale_t  =  11.07718404 
theta_graph  =  -1.0 
range_graph  =  0 
theta  =  0 

theta_l  =  acos  (cl  *  cos(theta)  /  cO) 

theta_2  =  acos  (c2  *  cos(theta)  /  cO) 

gO  =  ABS(  (cO  -  cl)  /  (dl  -  dO) ) 
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gi  =  ABS(  0.017  ) 

g2  =  ABS(  (cO  -  c2)  /  dO  ) 

range  =  c2/(g2*cos(theta_2))  *  (sin(theta_2)  -  sin(theta))  + 

+  cO/(gO*cos(theta))  *  (sin(theta_l)  -  sin(theta))  + 

+  h  •  cotan(theta_l)  + 

+  cl/(gl*cos(theta_l))  *  sin(theta_l) 

C  Vertical  dimension  for  each  graph  square  in  meters 
unit_dimension  =  range_scale  /  scale_r 

C  Dimensions  for  the  graph  for  the  full  profile  in  meters  and  radians 
range_14  =  2  *  range 

theta_14  =  theta_2 

DO  i  =  1,  100 

CALL  moveto_w(  theta_graph,  range_graph,  wxy  ) 
theta  =  theta +.0026179939 
theta_l  =  acos  (cl  *  cos(theta)  /  cO) 
theta_2  =  acos  (c2*cos(theta)/c0) 
range_graph  = 

+  c2/(g2*cos(theta_2))  *  (sin(theta__2)  -  sin(theta))  + 

+  c0/(g0*cos(theta))  *  (sin(theta_iy  -  sin(theta))  + 

+  h  *  cotan(theta_l)  + 

+  cl/(gl*cos(theta_l))  *  sin(theta_l) 

range_graph=(range_graph  -  range)/range_scale 
C  Scaling  factor  for  the  graph 

range_graph=  scale_r  *  range_graph 
theta_graph=  scalej  *  theta  - 1 
dummy  =  lineto_w(  theta_graph,  range_graph  ) 

END  DO 

Plot  the  non-mixed  layer  curve 

dummy  =  selcolor(  12  ) 
theta_graph  =  -1.0 
theta  =  0.0 
range_graph  =  0.0 
gO  =  ABS(  (cO  -  cl)  /  dl  ) 
theta_l  =  acos(crcos(theta)/c0) 

range  =  c0/(g0*cos(theta))  * 

+  (sin(theta_l)  -  sin(theta))  + 
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cO/(gl*cos(theta))  *  sin(theta_l  ) 


C  Dimensions  for  the  graph  for  the  bilinear  profile  in  meters  and  radians 
range_12  =  2  *  range 

theta_12  =  0 

DO  i  =  1,  100 

CALL  moveto_w(  theta_graph,  range_graph,  wxy  ) 
theta = theta + .0026179939 
theta_l  =  acos(cl*cos(theta)/cO) 
range_graph  =  c0/(g0*cos(theta))  * 

+  (sin(theta_l)  -  sin(theta))  + 

+  cl/(gl*cos(theta_l))  •  sin(theta_l ) 

range_graph  =  (range_graph  -  range)  /  range_scale 
C  Scaling  factor  for  graph 

range_graph  =  scale_r  *  range_graph 

theta_graph  =  scalej  *  theta  - 1 

dummy  =  lineto_w(  theta_graph,  range  graph  ) 

END  DO 
C 

C  Put  a  plus  for  the  estimated  minimum  range  and  its  angle 
C 

dummy  =  setcolor(  15  ) 
gamma  =  cl  /  cO 
G  =  1  +  gO  /  gl 

theta  =  SQRT(  (  1  -  gamma**2  )  /  (  G**2  -  1  ) ) 

range = (c0/(g0*cos(theta))) * ((1 + gO/gl ) *DSQRT( l-(cl/cO) *  *2  * 

+  cos(theta)**2)  -  sin(theta)) 

range_graph  =  ( range  /  range_12  )  -  1 
c  WRI1tE(*,*)  range_12,  range 

range_graph  =  scale_r  *  range_graph 
theta_graph  =  scale_t*  theta  - 1 
c  CALL  moveto_w(  theta_graph-0.05,  range_graph,  wxy  ) 
c  dummy  =  lineto_w(  theta_graph+0.05,  range_graph  ) 
c  CALL  moveto_w(  theta_graph,  range_graph-6^05,  wxy  ) 
c  dummy  =  lineto_w(  theta_graph,  range_graph+0.05  ) 

C 

C  Reference  information  for  the  graph 
C 

dummy  =  settextcolor(  13  ) 

CALL  settextposition(  2,  14,  curpos  ) 

WRITE  (string,  ’(F4.1)’)  unit_dimension  /  1852 
CALL  outtext(  ’Unit  dimensions:  ’) 
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CALL  settextpositipn(  3,  13,  curpos  ) 

CALL  outtext(  ’  1.5  deg  by’  //  string  //  ’nmi  ’) 
dummy  =  settextcolor(  12  ) 

CALL  settextposition(  4,  12,  curpos ) 

WRITE  (string,  ’(F5.i)’)  range_12  /  1852 
CALL  outtext(  ’Range  (0)  =  ’//  string//  ’nmi  ’) 

CALL  settextposition(  5,  12,  curpos  ) 

WRITE(  string,  ’(F4.2)’)  theta_12  *  57.29578 
CALL  outtext(  ’Theta  (0)  =  ’  //  string//  ’  ’ ) 

CALL  settextposition(  5,  32,  curpos  ) 

CALL  outtext(  ’degrees’) 
dummy  =  settextcolor(  14  ) 

CALL  settextposition(  6,  12,  curpos  ) 

WRITE  (string,  ’(F5.1)’)  range_14  /  1852 
CALL  outtext(  ’Range  (0)  =  ’/*/  string//  ’nmi  ’) 

CALL  settextposition(  7,  12,  curpos  ) 

WRITE(  string,  ’(F5.2)’)  theta_14  *  57.29578 
CALL  outtext(  ’Theta  (0)  =  ’  //  string//  ’  ’ ) 

CALL  settextposition(  7,  32,  curpos ) 

CALL  outtext(’degrees’) 
dummy  =  settextcolor(  15  ) 

CALL  settextposition(  8,  12,  curpos ) 

WRITE  (string,  ’(F5.1)’)  range  /  926 

CALL  outtext(  ’Range  (+)  =  ’//  string//  ’nmi  ’) 

CALL  settextposition(  9,  12,  curpos ) 

WRITE(  string,  ’(F5.2)’)  theta  »  57.29578 
CALL  outtext(  ’Theta  (  +  )  =  ’  //  string//  ’  ’ ) 

CALL  settextposition(  9,  32,  curpos  ) 

CALL  outtext(’degrees’) 
dummy  =  setcolor(  3  ) 

END 

CC 

CC 

CC  Subroutine  for  the  graph  description 

CC 

CC 

SUBROUTINE  graph_description(  cO,  dO,  cl,  dl,  d3,  h,  c2  ) 
INCLUDE  ’FGRAPH.FD’ 

REAL*8  cO,  dO,  cl,  dl,  c3,  d3,  h,  gl,  g2,  gO,  gain_bilinear, 

+  gain_mxdlyr,  theta,  RO,  Rl,  R2,  max_range,  CZ_width,theta_l, 
+  theta_2,  min_range,  temp,  c2 
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non  non 


CHARACTER’S  string 
RECORD  /videoconfig/  screen 
RECORD  /rccoord/  cuipos 
COMMON  screen 


Information  on  the  sound  speed  vs.  depth  graph 

dummy  =  settextcolor(15) 

CALL  settextposition(  2,  2,  curpos ) 

CALL  outtext(  ’  Depth, m  Sound  Speed,  m/s’) 

dummy  =  settextcolor(14) 

CALL  settextposition(  3,  2,  curpos  ) 

WRITE  (string,  ’(I5)’)INT(  dO  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  3,  15,  curpos ) 

WRITE  (string,  ’(I5)’)INT(  cO  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  4,  2,  curpos ) 

WRITE  (string,  ’(I5)’)INT(  dl  ) 

CALL  outtext(  ”  //  string  //  ’  ’) 

CALL  settextposition(  4, 15,  curpos ) 

WRITE  (string,  ’(I5)’)INT(  cl  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  5,  2,  curpos  ) 

WRITE  (string,  ’(I5)’)INT(  dl  +  h  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  5,  15,  curpos ) 

WRITE  (string,  ’(I5)’)INT(  cl  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  6,  2,  curpos ) 

WRITE  (string,  ’(I5)’)INT(  d3  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  6,  15,  curpos ) 
c3  =  cl  +  0.017*(d3-dl-h) 

WRITE  (string,  ’(I5)’)INT(  c3  ) 

CALL  outtext(  ’  ’  //  string  //  ’  ’) 

Gain  information,  bilinear  then  mixed  layer 

temp  =  temp 

gain_mxdlyr  =  0 

gO  =  ABS(  (cO  -  cl)  /  (dl  -  dO)  ) 
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gl  =  ABS(  0.017  ) 
g2  =  ABS(  (cl  -  c2)  /  dO  ) 


gain  bilinear  =  10  *  LOG10(  (2  *  SQRT(  (l+g0/gl)**2-l  ))  / 

+  ■  ( (1+gO/gl)  -  SQRT(  (l+g0/gl)^*2-l  ))) 

dummy  =  settextcolor(12) 

CALL  settextposition(  8,  2,  curpos ) 

CALL  outtext(  ’  CZ  Gain  =’) 

CALL  settextposition(  8,  12,  curpos ) 

WRITE  (string,  ’(F'^-l)’)  gain_bilinear 
CALL  outtext(  ’  ’  //  string  //  ’  ’) 

CALL  settextposition(  8,  18,  curpos ) 

CALL  outtext(  ’db’) 

gain  mxdlyr  =  10  *  LOG10(  (2  *  SQRT(  (l+g0/gl+g0/g2)**2-l  ))  / 
+  "  ( (l+g0/gl+g0/g2)  -  SQRT(  (l+g0/gl+g0/g2)**2-l  ))) 

dummy  =  settextcolor(14) 

CALL  settextposition(  10,  2,  curpos ) 

CALL  outtext(  ’  CZ  Gain  =’) 

CALL  settextposition(  10,  12,  curpos ) 

WRITE  (string,  ’(F4.1)’)  gain_bilinear 
CALL  outtext(  ’  ’  //  string  //’’  ’) 

CALL  settextposition(  10,  18,  curpos ) 

CALL  outtext(  ’db’) 

END 
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APPENDIX  F 


USER’S  MANUAL  FOR  CZ  GAIN 


CZGAIN  requires  an  MS-DOS  based  system  with  high  resolution  graphics  capability 
(VGA  or  better).  The  model  was  developed  using  ray  tracing  rather  than  normal 
mode  theory. 

The  user  inputs  the  following  data: 

•  Surface  sound  speed,  m/s 

•  Depth  of  mixed  layer,  m 

•  Sound  speed  at  bottom  of  mixed  layer,  m/s 

•  Minimum  sound  speed  below  the  mbced  layer,  m/s 

•  Depth  of  that  sound  speed,  m/s 

•  Bottom  depth,  m 

While  the  bottom  sound  speed  gradient  is  fixed  at  0.017  /s,  the  bottom  depth  is  used 
when  drawing  SSP’s. 

The  surface  sound  .speed  and  minimum  sound  speed  will  always  be  less  than  the 
sound  speed  at  the  bottom  of  the  mixed  layer. 

The  depth  of  the  mixed  layer  will  be  no  less  than  1  meter. 

The  thickness  of  the  isospeed  layer  will  be  zero  or  greater. 

When  developing  the  keys  to  use  for  the  data  manipulation,  efforts  were  made  to  tie 
each  data  point  to  a  specific  mnemonic  (e.g.  Thickness  of  isospeed  layer  changed  by 
t).  No  standard  convention  could  be  established.  The  numbers  work  from  left  to 
right  on  sound  speed  and  depth  pairs.  Where  no  depth  or  sound  speed  was  used 
(e.g.  surface  depth  is  always  zero)  that  key  was  used  for  the  next  parameter.  The 
numbers  always  decrease  the  respective  value  and  using  the  SHIFT  key  with  the 
number  always  increases  the  respective  value.  The  letters  were  added  at  the 
programmers  personal  convenience.  No  enter  key  is  required,  the  program 
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automatically  captures  the  keystroke  and  recomputes  the  output.  The  following  is 
a  list  of  keys  which  will  manipulate  the  data. 

•  1  -  Decreases  surface  sound  speed 

•  !,  A,  a  -  Increases  surface  sound  speed 

•  2  -  Decreases  sound  speed  at  bottom  of  mixed  layer 

•  S,  s  -  Increases  sound  speed  at  bottom  of  mbced  layer 

•  3  -  Decreases  depth  of  mixed  layer 

•  #,  D,  d  -  Increases  depth  of  mbced  layer 

•  4  -  Decreases  minimum  sound  speed 

•  $,  F,  f  -  Increases  minimum  sound  speed 

•  5  -  Decreases  depth  of  minimum  sound  speed 

•  %,  G,  g~  Increases  depth  of  minimum  sound  speed 

•  6  -  Deceases  thickness  of  isospeed  layer 

•  H,  h  -  Increases  thickness  of  isospeed  layer 

•  7  -  Decreases  bottom  depth 

•  &,  J,  j  -  Increases  bottom  depth 

•  Q,  q,  E,  e  -  Exits  the  program 

If  the  program  is  modified,  the  Key  Entry  Functions  (if  using  Microsoft  Fortran) 
patch  will  need  to  be  installed  into  the  proper  libraries  before  compiling. 
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